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Abstract: We investigate the relation between the annihilation of dark matter (DM) 
particles into lepton pairs and into 2-body final states including one or two photons. We 
parametrize the DM interactions with leptons in terms of contact interactions, and cal¬ 
culate the loop-level annihilation into monochromatic gamma rays, specifically computing 
the ratio of the DM annihilation cross sections into two gamma rays versus lepton pairs. 
While the loop-level processes are generically suppressed in comparison with the tree-level 
annihilation into leptons, we find that some choices for the mediator spin and coupling 
structure lead to large branching fractions into gamma-ray lines. This result has implica¬ 
tions for a dark matter contribution to the AMS-02 positron excess. We also explore the 
possibility of mediators which are charged under a dark symmetry and find that, for these 
loop-level processes, an effective field theory description is accurate for DM masses up to 
about half the mediator mass. 
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1 Introduction 

While there is very strong evidence for the existence of dark matter (DM) on the scales of 
galaxies, galaxy clusters and cosmology [1, 2], determining its specific particle properties 
remains one of the most important open questions in modern physics [3]. There are many 
elementary particle candidates for the dark matter in the universe (see e.g. [4] for a review 
of particle DM candidates), but one of the most compelling class of models is based on a 
striking coincidence: a stable particle with approximately weak-scale mass and couplings (a 
so-called weakly-interacting massive particle, or WIMP) will, assuming a standard thermal 
history of the universe, possess a thermal relic density today that roughly corresponds to 
that of DM. This coincidence is known as “WIMP miracle”. 

Much effort is under way in searching for WIMPs. The most direct probe of WIMPs as 
DM candidates is the search for scattering of Galactic WIMPs with ordinary matter. Many 
direct detection experiments have searched for this signal, with no dehnitive detections as 
of yet [5, 6]. The bounds from such searches have begun to have significant implications for 
our understanding of viable WIMP DM candidates, ruling out the most naive realizations 
of the WIMP paradigm. Of course, not all DM interactions can be effectively probed by 
these experiments. 

A second technique to search for WIMPs is to attempt to produce them in high-energy 
particle collisions. Once produced, a WIMP would escape undetected from the experiment, 
leading to a missing energy signature [7]. Searches for missing energy have a very long 
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history, and have been designed for a number of purposes, besides models featuring a stable 
WIMP. Recently, a focus on comparing the bounds on interactions of WIMPs with other 
particles from these searches to indirect and direct detection experiments has emerged. This 
exercise started with a parametrization of DM interactions by contact interactions [8-10], 
and has evolved to include various simplified models which feature a collider-accessible 
mediator particle as well [11]. Many collider analyses are now being designed explicitly 
with the purpose of searching for DM [12]. 

A third technique, known as indirect detection, searches for the visible particles that 
result from annihilations of WIMPs in the Galaxy or in outer systems [5] . These annihila¬ 
tions are rare enough that they do not lead to any significant depletion of the DM density, 
but they can potentially produce measurable signals. Many experiments have sensitivities 
good enough to probe annihilation rates that correspond to the rate needed in the early 
universe in order for the WIMP Miracle to occur [13]. We must keep in mind, however, 
that the kinematics of annihilation at the present time are very different from those needed 
in the early universe, and thus the annihilation rate may be significantly different as well 
[14]. 

A variety of target Standard Model (SM) particles can be employed for the purposes 
of indirect searches. Neutrino telescopes have been used to search for the annihilation of 
WIMPs captured in the Sun, photons from radio through gamma rays have been used to 
place bounds on DM annihilation, and cosmic rays can also be employed to search for a 
WIMP signal. The most promising cosmic rays for such purpose are antimatter particles, 
as we believe that ordinary astrophysical processes do not produce antimatter in great 
abundances, and that therefore the only background should be secondary, i.e. arising from 
cosmic ray interactions with interstellar material. Such secondaries are generally much less 
abundant and have much less energy than primary cosmic rays, so that WIMP annihilation 
could potentially yield a detectable primary contribution. 

In 2008, the PAMELA experiment precisely measured the positron fraction in the 1.5- 
100 GeV range [15]. The positron fraction was found to sharply increase around 5 GeV in a 
manner completely inconsistent with secondary sources in essentially all predictive cosmic- 
ray models [16, 17]. The Fermi Large Area Telescope (LAT) later confirmed that the 
positron faction continues to rise in the 100-200 GeV range [18]. The International Space 
Station-based experiment Alpha Magnetic Spectrometer (AMS-02) extended measurements 
to 350 GeV [19]. Their 6.8 x 10® electron and positron events confirm the PAMELA excess 
with high precision. The slope of the AMS-02 spectrum decreases by an order of magnitude 
between 20 and 250 GeV and has no measurable fine structure or anisotropy in the arrival 
direction. The second AMS-02 data release last fall extended the analysis to 500 GeV, and 
found that above ~ 275 GeV the positron fraction ceases to increase [20]. 

These experiments imply that there must be another significant source of high energy 
positrons besides secondary processes. While a Galactic WIMP annihilation origin is still 
possible [21], nearby pulsars have been proposed as another possible explanation for the 
positron excess [22, 23]. 

WIMP annihilation into charged states also generically produces photons via loop-level 
processes or from final-state radiation. In the latter case, the secondary photons have an 
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Figure 1. Feynman diagrams for the processes yy —?► // and yy —>■ 77 . 


energy spectrum peaked towards low energies. In the former case, however, when DM 
annihilates at close to zero relative velocity into 7 X, conservation of energy gives a photon 
with a fixed (monochromatic) energy of 

=“x (‘ - £|) ■ 

where is the WIMP mass. A gamma-ray line from these processes could stand out from 
the continuous astrophysical background for high enough values of and be detected by 
HESS or Fermi LAT. Since gamma rays propagate in straight lines from their production 
location, the line of sight of their spatial origin can easily be determined and compared 
with expectations based on various WIMP models and DM density prohles. While there 
was recent excitement about a 130 GeV line in the public Fermi LAT data [24], the signal 
is not significantly present in the Fermi Collaboration’s latest analyses based on their Pass 
8 software upgrade, and no other spectral lines have been detected [25]. 

In this paper we study the monochromatic gamma-ray yield stemming from DM models 
where the DM particle annihilates into a charged lepton-antilepton pair. In particular, we 
carry out a model-independent study based on contact interactions in the context of an 
effective field theory (EFT) description of the DM-lepton interaction vertex. For a given 
contact interaction type, we essentially simply “close the fermion loop” as shown in figure 1 
and calculate the cross section for a given annihilation cross section into leptons 

As a result, we calculate in a model-independent way the monochromatic 
gamma-ray yield for e.g. theories that purportedly explain the measured positron excess; 
vice versa, our results allow to compute (cru)^^_^jj from a given measured gamma-ray line 
and model to explain it. 

The remainder of this paper is organized as follows: In section 2, we describe the EFTs 
structure under consideration. We calculate the ratio in section 

3.1. In section 3.2, we comment on the applicability of the EFT assumption in the case 
of mediators charged under the dark symmetry. Before concluding, as an example of the 
utility of our ratio results, we compute the rate of annihilation into monochromatic gamma 
rays in section 3.3, assuming a WIMP explanation for the positron excess. 

2 Model Setup 

We consider a model of WIMP DM where the WIMP is a SM gauge singlet Dirac fermion 
X, and we assume that the new particle mediating the interactions between the WIMP and 
SM helds is heavy enough to be “integrated out”. The natural way of parametrizing these 
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Name 

Operator 

X 

VV 

AV 

VA 

X7^7®x/7m/ 

X7^x/7p7V 

X 

AA 


7 

RV 

TV 

Xl^PRX^liJi^ 

Xl^Plxflp^f 

z 

RA 

LA 

xx^Pnxht^x^f 

7 

VR 

VL 

xl^xh^lPRf 

Xl^xfli^Pbf 

z 

AR 

AL 

xx^x^xhi^PRf 

Xl^l^xfltiPbl 

7 


Name 

Operator 

X 

SS 

xxff 


PS 

xi^xff 


SP 

xxfi^f 

7 

pp 

Xl^xfl^f 


RR 

Xl^PRXhiiPRf 


RL 

Xl^PRXfli^Pbf 


LR 

Xl^PLXfli^PRf 

7 

LL 

xYPlxflfiPlf 



Table 1. Effective operators for WIMP-SM interactions considered in this work. The X column 
indicates the second gauge boson in the final state of the XX process. 


interactions, then, is to have non-renormalizable contact operators suppressed by a mass 
scale that is related to the mediator’s mass and couplings. In this work we focus on interac¬ 
tions between DM and a SM fermion /. The lowest-dimension, and thus least suppressed, 
operators are therefore dimension six products of DM bilinears and SM bilinears. While 
we could consider operators containing mixed bilinears as well, these are all equivalent 
through Fierz transformations to linear combinations of the class we do consider. 

The interactions are chosen to be of the form [Trx][/S/], where T and S are matrices 
in the set { 1 , 7 ^, 7 ^, 7 ^ 7 ^, 7 ^Pr, with Pr = ^(1 + 7 ®) and Pl = ^(1 — 7 ®) being 
the usual helicity projectors. This is obviously an over-complete basis of operators, as the 
vector and axial interactions can be constructed out of the chiral interactions or vice versa. 
However, due to the chiral nature of the SM gauge interactions the chiral operators may 
well be of interest, despite their linear dependence on the others. The full list of EFT 
interactions we consider are shown in Table 1. Our naming scheme is simple. For scalar 
operators, S and P stand for scalar (1) and pseudoscalar ( 7 ^) couplings, respectively. V, 
A, R and L identify vector ( 7 ^), axial ( 7 ^ 7 ^), right-handed {x^Pr) and left-handed {x^Pl) 
couplings. 

The quantity of interest here is the cross section ratio 


Rxj 
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where X is another SM vector boson. ^ As discussed before, and are 

computed by evaluating the diagrams shown in figure 1. We take A = 7 for the scalar-type 
EFTs to maximize photon production. For operators of the form [x^^x\[f'P‘^lf]^ X = Z 
unless the WIMP and SM pieces both contain a factor of 7 ®, in which case we also take 

^Note that the cross sections in this work are calculated at a specific value of the WIMP relative velocity 
V, as opposed to being thermally averaged. 
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X = 7- See appendix B for details behind these choices, which are related to Yang’s 
theorem. 

It is important to note that, while we use the language of EFT interactions, any s- 
channel mediator with appropriate couplings will in fact give an identical cross section ratio 
to our EFT-based approach, as the s-channel propagator in both diagrams is identical and 
the mediator cannot carry the loop momentum. Thus, our resulting ratios are applicable 
for these theories beyond the regime where the EFT has captured all DM physics of interest. 
After presenting our EFT results we will discuss the case of t-channel mediators in some 
detail, determining where the EFT analysis holds for Rxj and where it breaks down and 
the full set of box diagrams must be calculated explicitly. 

The AA, AR, AL and AV operators have p-wave suppressions at tree level. This implies 
that the annihilation rate in the Milky Way halo today is ~ 10““^ that in the early universe, 
at the time when a thermal WIMP would freeze out. Since the cross section required today 
to explain the positron excess is instead ~ 10^ times the thermal relic one, this potential 
explanation of the positron excess requires a non-thermal cosmology, e.g. DM production 
in moduli decays [26] or a modihed expansion history, e.g. due to quintessence [27]. 

These suppressions also mean that we cannot consistently take the relative WIMP 
velocity to be zero when computing Rxj for all of our theories, as is often done. We 
found that the 3-body processes XX ffl and XX ff^ did not significantly alleviate 
the p-wave suppressions, nor the additional helicity suppression found for in the 

A A theory, demonstrating that gamma ray lines are the natural signals for our EFTs. 

We take the velocity for the tree-level process to be r'x,r/c = 10“^, a choice motivated 
by the typical halo velocity dispersion [28]. While gamma-ray line searches focus on the 
Galactic Center, the velocity dispersion in this region is within a factor of two of v-^^r [29], 
so we are justified in using the standard approximation = '^x calculate 

((Tu)xx^> 77 - Accounting for the effect of Sgr A* on the WIMP velocity dispersion is beyond 
the scope of this paper. 

3 Results 

3.1 Cross Section Ratios 

We calculate the ratio Rx,t in each of the EFTs in Table 1 using the Mathematica pack¬ 
ages FeynArts [30], FormCalc and LoopTools [31]. The spin-averaged cross sections are 
computed with WIMP velocity v^jc = 10“^ in the center of mass frame. We collect our 
analytic results for the relevant tree-level cross sections in appendix A.l. The tree-level and 
loop-level analytic amplitudes for the SP and AA interactions found in appendix A.2 and 
A.3 were found to agree with the numerical results. Notice that while our EFT interactions 
are non-renormalizable, the one-loop amplitudes are finite. 

The cross section ratios i?x,r for a representative set of our EFTs are shown in figure 
2 with WIMP masses ranging from ~ 100 GeV to ~ 10 TeV. Ratio curves are omitted for 
some EFTs as they are related simply to others shown in the plot: 

• In the case of the scalar theories, the WIMP and fermion (gamma) parts of the tree- 
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Selected cross section ratios 



(GeV) 


- AA (77) 

- AR (77) 

- AV (7Z) 

- RA (77) 

RR (77) 

- RV (jZ) 

- VV (7Z) 

- SS (77) 

'— VR (7Z) 

VA (7Z) 


Figure 2. Cross section ratios for the representative set of the EFTs discussed in section 3. The 
particle name following each EFT label indicates the choice of X that maximizes photon production 
through the loop diagram. For theories with X = Z, curves begin at the Z production threshold 

= Mzl2 = 45.6 GeV. 

level (loop-level) cross sections decouple. Since the WIMP components are the same 
for the tree and loop-level cross sections, the SS and PS results are the same, as are 
the PS and PP ones. Moreover, the loop parts of the amplitudes differ by a factor of 
order 1 for scalar and pseudoscalar coupling to /. We therefore omit the PS, SP and 
PP results, as they are numerically very similar to the SS one. 

• Averaging over initial WIMP spins allows us to omit LA, LV, LL and LR from our 
plot, since they give the same results as the RA, RV, RL and RR theories. 

• The AR, AL, RR and RL theories allow annihilations of WIMPs to two photons. 
Since photons couplings to / are helicity independent, we also left AL and RL out of 
our plot. 

• The Rz,t ratios for the VR and VL theories differ by a factor ^ J^+sfiP e ) ~ 
since left and right-handed rs couple differently to the Z, so we only show the former. 

Most of the EFTs give Rx,t ~ 10“^° — 10“®. However, for AA, AL, AR, Rx,t ~ 
10“^ — 10 and for AV, Rz,t ~ 1- While surprising at first that a loop-level to tree- 
level diagram give such close cross sections, this is easily understood by examining the 
relevant suppressions. As mentioned before, the tree-level cross section for the AV theory 
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Figure 3. The box diagram contributing to (av)^^-yjj in our theory containing a t-channel 
mediator r. 

shown in appendix A.l is p-wave suppressed, but the loop one is not, leading to a large 
ratio. Similarly, while the AA, AL and AR theories are p-wave and helicity suppressed 
at tree-level, these suppressions are alleviated in the loop diagrams. The ratios for these 
three theories demonstrate the transition between the regime in which the tree-level cross 
section is dominated by the helicity suppressed term (M^ ^ 3 x 10^ GeV) and the p-wave 
suppressed term (M^ ^ 3 x 10^ GeV). In the case / = // or e“, the tree-level helicity 
suppressed terms becomes negligible for these theories, and their ratios are constants equal 
to their asymptotic values of 1 — 10. 

3.2 EFT for t-channel Mediators 

The loop diagrams we consider require care as arbitrarily large loop momenta run through 
our effective operators. As mentioned before, if the contact operator comes from an s- 
channel mediator, the mediator is not contained in the loop and the EFT gives the correct 
result for the ratio of cross sections. In the case of a charged t-channel mediator, we must 
consider the box diagram shown in figure 3 (along with the crossed ones). In this case, the 
mediator’s propagator could lead to significant corrections to the cross section from the 
EFT prediction at high loop momenta. 

To explore this potential issue, we considered a theory in which the WIMP-SM inter¬ 
action is of the form 


Ant = [x(As + Ap 7 ^)r]f+ -b {h.c.), (3.1) 

where y is the WIMP (taken to be Majorana for simplicity) with mass rriy^ and f~^ is a 
charged scalar with mass . This is a simplified model of a SUSY theory in which one tan 
slepton is the dominant contributor to DM annihilations and the DM neutralino’s gauge 
eigenstate content is unknown. We consider the CP-conserving case, which requires that 
Xs, Ap G M. 

To lowest order, the effective Lagrangian for this theory will contain a linear combina¬ 
tion of four-fermion operators from Table 1. These are found by applying Fierz transfor¬ 
mations to the matrix element for xx , giving 

Ant ^ Aflf = {[xx][fT] + [xi^x][fl^r]) - [t7m7^t-] 

- ^^^[X7Vx][r7/.r]- (3.2) 
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Figure 4. and plotted as a function of for As = 1 and Ap = 0. The 

agreement between the two cross sections is better for any other choice of the coupling constants. 

Note that in the case of a Dirac WIMP, there would be additional terms involving its tensor 
and vector currents. In the chiral case |As| = |Ap|, only the AA and AV terms are present, 
implying that will be large for this theory. In what follows we will call the EFT and 
exact cross sections and The latter has been calculated in many 

papers (see e.g. [32-36]) in the = 0 limit. Here we use the result from [35] which is 

corrected in the footnote on page 8 of [36]. 

In figure 4 we plot and as a function of m-^ for ttit = 10 TeV, 

A 5 = 0 and Ap = 1. The cross sections agree to better than 1% all the way up to 

~ ^rrif. Near ~ ruf the EFT underestimates the cross section when the slepton 
can be put on shell. For my_ > mf the EFT result is larger than the exact result since the 
EFT loop only contains three propagators. The coupling constants’ values have no effect 
on this agreement. It is surprising that the EFT is valid for such high values of and 
furthermore that it works so well for our SUSY-like theory, as this is the classic internal 
bremsstrahlung case in which EFT is usually expected to be inaccurate. 

Finally, we consider whether our calculation of changes in this t-channel theory. 
The cross section gets a new contribution from the internal bremsstrahlung 

diagram for XX ffl- This diagram only contributes significantly to if the 

helicity suppression in the lowest-order process xX ff is very large. Since we are 
interested in the final state f = t, this is not the case [33], and we can safely ignore 
the internal bremsstrahlung diagram for the range of WIMP masses considered above; 
for the largest WIMP masses considered the helicity suppression is comparable to the 
one-loop suppression, leading to a change of only a factor of a few in the ratio. Our 
type of EFT analysis of Ry^r from above thus applies for , although of course 

having to consider a linear combination of effective operators in eq. (3.2) adds additional 
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Figure 5. (M^, contours for the AR theory predicted using the calculated ratio 

and preferred regions at 99% CL obtained from positron data in [38], for Vy^lc= 10“^. 

complexity. For a lighter hnal state fermion or much heavier WIMP, a more careful analysis 
of including the internal bremsstrahlung diagram would have to be performed, 

but note that the WIMP mass for a simplified theory of this type has been shown to have 
an upper bound comparable to the masses considered here [37]. 

3.3 Predicting the Gamma Ray Flux 

As an example of the utility of our results, suppose now we postulate a WIMP annihila¬ 
tion explanation for the rising positron fraction. Analyses of the hrst release of AMS-02 
data found that TeV-scale WIMPs annihilating to rr provide the best explanation of the 
positron excess in terms of DM annihilating to SM particles [38-40] . This has not changed 
with AMS-02’s latest data release [41]. To make quantitative predictions for the gamma- 
ray line yield, we use for reference the allowed regions in the (m^, plane from 

[38]. Our calculated ratios enable us to predict the gamma-ray monochromatic line signal 
induced through the loop diagram for a given contact interaction. 

By rescaling the 99% CL allowed regions in the (M^^, plane given in [38] by 

Rx,t and putting in a factor of 1/2 for theories with X = Z, we calculated the regions in 
the (M^, plane preferred by a WIMP annihilation explanation to the positron 

fraction anomaly. Figure 5 shows the resulting regions for the AR EFT along with recent 
upper limits on from HESS [42] and Eermi LAT [43] line searches. The pre¬ 

dicted allowed regions for AA, AL and AV he above those for AR while the regions for 
all other theories fall well below current observational constraints, by at least 3 orders of 
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Figure 6. Same as in figure 5, but with increased (decreased) by an order of magnitude in the 
left (right) plot. 

magnitude. HESS-II, GAMMA-400 and CTA line searches are only expected to tighten 
these constraints by an order of magnitude [44] , and so will not be able to probe line signals 
from AMS-fitting points in these theories. 

For comparison, we also show the latest constraints on {av)-^^-^rT based on CMB data 
from Planck and other experiments [45] as well as Fermi’s observations of diffuse gamma- 
ray emission in the Milky Way [46], which are in tension with WIMP explanations of 
the positron fraction. Note that the rescaled CMB constraints are quite conservative at 
high M^. In this part of parameter space, the p-wave term dominates the tree-level cross 
section, which means the limits are a factor 1 larger, where Ux,recom is the 

characteristic relative WIMP velocity at recombination. 

Figure 6 shows the velocity dependence of the results for the AR theory. The differences 
between the two panels stem from changes in the relative strengths of the helicity and p- 
wave suppressed parts of the tree-level cross section. While increasing or decreasing by 
an order of magnitude changes the shape of the limit regions, it does not affect whether 
they fall above or below the line search constraints for any of our theories. We therefore 
find that AA, AR, AL and AV contact interactions cannot explain the positron excess as 
they would violate bounds from gamma-ray line searches, in addition to being in tension 
with Fermi diffuse and CMB constraints. All other contact interactions remain viable. 

4 Discussion and Conclusions 

In this study we have calculated the cross section for Dirac WIMPs annihilating to monochro¬ 
matic gamma rays stemming from an annihilation mode to a charged lepton pair. To relate 
and we used an over-complete set of dimension 6 operators. The ra- 
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tio Rx,t varies by many orders of magnitude depending on the theory under consideration, 
since helicity and p-wave suppressions in the tree-level cross section are sometimes allevi¬ 
ated at one loop. We also investigated the validity of the EFT description for a specific 
SUSY-inspired case (t-channel stau-mediated annihilation into tau pair) and demonstrated 
that the EFT accurately describes theories where a charged boson mediates SM-WIMP 
interactions. Since the EFT gives the same value for Rx,t as an s-channel mediator, this 
shows that our analyses hold for a wide range of theories whose UV completion involves a 
mediator heavier than the WIMP. 

We used our cross section ratios to predict the gamma ray signal given an annihilating 
WIMP explanation for the positron excess. Describing the positron excess using contact 
operators involving an axial WIMP bilinear (namely the AA, AR, AL and AV operators) 
gives a monochromatic gamma-ray line orders of magnitude brighter than current obser¬ 
vational limits from Fermi LAT and HESS. Our approach is complementary to constraints 
on diffuse gamma ray emission and CMB measurements, and will prove useful in directing 
gamma ray searches should a signal from WIMPs annihilating to leptons be detected in 
the future. 
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A Analytic results 

A.l A selection of tree-level cross sections 

Table 2 shows the tree-level cross sections for several theories, evaluated to order v^. 


Theory 

{(jv) expanded to 0{v‘^) 

SS 

X^vl 2 2J 

PP 


AV 


AA 

\haj ( 1 

47r y Ml) ^ ^ ^ + ^2 J 


Table 2. Tree-level cross sections for a representative subset of EFTs. 


A.2 Loop cross section for SP theory 

The total loop amplitude is found by summing the diagrams shown in figure 7: 


-^loop = A[n-/(fe 2 )u^(A:i)]e* (fe3)<(A:4)iV'^"(A:3, ^ 4 ), (AT) 

where A is the coupling at the blob and A^'^(/c 3 , k^) = N^'^{k 3 , ^ 4 ) -|-iV^^(/c 4 , k^) is the sum 
of loop integrals for the crossed and uncrossed diagrams. The former is given by 


N^‘'{k3,k4) = i-ie)^ 


tr[(yi - 1^4 -b mf)Y{p + + h+ ^fh^] 

(27r)"^ [{p — A:4 — mj] — mj] [{p + k^Y — m^] 


(A. 2 ) 


The trace evaluates to Airufe^’^'^^ksctk^p by standard identities. After introducing Feynman 
parameters, the integral becomes 


N^''{k 3 ,kY = j dx j dz 


dH 


— ruj + 2xzk3 ■ kYp 


I = p — xk^ + zk3. 

(A.3) 
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Wick rotating and integrating yields 


N^''ik3,h) 




3aM/3 


Li 2 




(A.4) 

This can be converted to the more pleasant expression found in [CITE] by applying Lan- 
den’s identity: 




log^ 



s + 


yjs{s — Amj) 
2m‘j 



(A.5) 


Since this expression is symmetric under (/i ^ ^ 4 ), the total amplitude is found to 

be 


■^loop = w{k2)u\ki)]e^^^‘^^£*^{k3)el{ki)k3o,kii3 log^ 


s + s{s — 4m‘j) 


2m‘j 


The spin and polarization-averaged matrix element is therefore 

s + ^Js{s — 4mj) 




2m?j 


(A. 6 ) 


(A.7) 


Doing the Lorentz invariant phase space integration gives the final result for the cross 
section: 

2 

^ ^ 2 / jO\ i i-»l ! r-i { /-> A Tv-fco 2 ^ \ 

"a ml 


{av) 


SP 




1 — 

1 

log^ ( 

\ ® / 

\ 


xx^ii 2567r4 

A.3 Loop cross section for AA theory 

Using the result for M.z^'yy from [CITE], we find 


s + 


2m?j 


(A. 8 ) 


{av) 


XX^ll 


X^a^M‘^s^ 

dvr^ 


m 


IH- - log^ 

s 


s + s{s — Am'j) 


2mj 


(A.9) 


Note that equations 3.24 and 3.25 in [CITE] are both missing a factor of M 


-1 


B Yang’s theorem 

Yang’s theorem [47] states that a massive spin-1 particle cannot decay into a pair of identical 
massless spin -1 particles when the particles are all on shell, which has implications about 
{av)-^-^-^r^^ for vector-type EFTs. Any of our effective operators can be interpreted as the 
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limit of a Z' theory as M^/ —)• cx). Here we will use this perspective to make some remarks 
about the process xx —)• 77 for vector-like theories of the form [xr^x][^ 5 ]^r]. 

For such a theory, the EFT amplitude for this process can be written schematically as 

= [vT^u] • (B.l) 

where N^p^j represents the loop piece of This is just the limit of the following 

amplitude in a Z' theory: 

~ ■ ^2 _ Xf2 ~ ]^2 ^ ■ ^vpa ■ Ep^o-) (B.2) 

where q is the sum of the WIMP momenta. When the Z' is on shell, the factor in 
the propagator with Lorentz indices is the projection onto the its physical polarizations, 
— Yle^‘qfj,=o Thus Yang’s theorem requires that {q^ — vanish when 

q = y/s = Mz'. 

There are two ways this can happen. When = 0, this condition will 

of course be satisfied. This is true when = 7^, where Furry’s theorem gives the much 
stronger result that Wpo- and thus = 0 for all q. \M.^-^-^^^{q)\^ also vanishes 

for the VA theory. The only other way this can happen is if contains a factor 

of to cancel the one in the propagator. This only occurs when the Z' has nonzero 

axial coupling to the WIMP and fermion. We have verified this using FeynCalc and the 
analytic expression for Z' —>■ 77 [ 48 ] . 

However, [ 48 ] demonstrated that the triangle diagram for Z' —)• 77 gives rise to an 
anomaly, which is cancelled by mixing between the Goldstone bosons arising from sponta¬ 
neous symmetry breaking and the Z'. This means that vanishes for Z' models? 

Of course, effective operators with axial coupling to the WIMP and fermion can arise in 
theories such as the t-channel one we considered in section 3 . 2 . In this case there is no 
anomaly since no gauge symmetry is being violated, and is nonvanishing. 


^Since the matrix element in this cross section is zero, the third case discussed in [49] cannot occur. 
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